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ABSTRACT 

Aims. In this pilot study, we examine molecular jets from the embedded Class I sources, HH26 and HH72, to search, for the first 
time, for kinematic signatures of jet rotation from young embedded sources. 

Methods. High-resolution long-slit spectroscopy of the H2 1-0 S(l) transition was obtained using VLT/ISAAC. The slit was placed 
perpendicular to the flow direction about 2" from the sources. Position- velocity (PV) diagrams are constructed and intensity- 
weighted radial velocities transverse to the jet flow are measured. 

Results. Mean intensity- weighted velocities vary between v LSR ~ -90 and -65 km s -1 for HH26, and -60 and -10 km s _1 for HH 72; 
maxima occur close to the intensity peak and decrease toward the jet borders. Velocity dispersions are ~ 45 and ~ 80km s for 
HH 26 and HH 72, respectively, with gas motions as fast as - 100 km s -1 present. Asymmetric PV diagrams are seen for both objects, 
which a simple empirical model of a cylindrical jet section shows could in principle be reproduced by jet rotation alone. Assuming 
magneto-centrifugal launching, the observed HH 26 flow may originate at a disk radius of 2-4 AU from the star with the toroidal 
component of the magnetic field dominant at the observed location, in agreement with magnetic collimation models. We estimate 
that the kinetic angular momentum transported by the HH 26 jet is ~ 2 x 10~ 5 M Q yr~' AU km s . This value (a lower limit to the 
total angular momentum transported by the flow) already amounts to 70% of the angular momentum that has to be extracted from 
the disk for the accretion to proceed at the observed rate. 

Conclusions. These results of this pilot study suggest that jet rotation may also be present at early evolutionary phases and support 
the hypothesis that they carry away excess angular momentum, thus allowing the central protostar to increase its mass. 

Key words, star formation - jets and outflows - angular momentum transport 



. ■ 1. Introduction cess angular momentum from the circumstellar disk. This is 

required to permit material to fall to low angular momentum 

O ; The appearance of jets from young stars and active galac- orbks and thus be finally accreted ontQ the protostari thereby 

OO , tic nuclei is arguably one of the most enigmatic and strik- increasing its mass up t0 its final value In addition , jets have 

O ■ ingly beautiful phenomena in astrophysics. Yet we understand the important role of dispersing the Mailing material from the 

> ; very little about the detailed physical processes which con- surrounding circumst ellar envelope and hence influence the fi- 

. rH , tnve to produce these jets and still debate the general mech- nal mass of the central protostar 
■ anisms responsible for them (for the mos t recen t reviews see 

Jh I iBally. Reipurth & DavisL l2Q07t iRay et all I2007I) . The inflow Recent observations give evidence of a small difference 

. . , and outflow of material is believed to be mediated by the in- between the radial velocities of lines emitted from two op- 

teraction of magnetic and centrifugal forces in the close envi- posed sides of a jet (with respect to the symmetry axis). This 

ronment of a newborn star. In particular, a magnetised accre- was interpreted in the context of gas rotation. The first ten- 

tion disk channels material towards the protostar (on which it tati ve results came fro m observations of the HH212 outflow 

will finally be accreted via coronal loops) and provides the nec- by iDavis et al.l (l2000l) . Velocity shifts of a few km s~' were 

essary magnetic forces to launch material off the disk surface found for the H2 lines across the breadth of the jet. However, 

into a jet/outflow. One fundamental problem to which jets are the rather coarse spatial and spectral resolution of the data pre- 

believed to provide a solution, is the necessary removal of ex- vented any definitive conclusions, especially given that the re- 

gion of the jet investigated in that study is located at not less 

Send offprint requests to: A. Chrysostomou than 10 4 AU from the source, too far from the acceleration zone 



and then ICoffev et alj (12004 120071) and IWoitas et al.1 ([2005) 
presented better evidence for jet rotation from several small- 
scale jets emanating from classical T-Tauri stars, including DG 
Tau, RW Aur, Th 28 and CW Tau. These high angular resolu- 
tion optical and near-ultraviolet observations, with slits placed 
parallel and perpendicular to the flow axis, are from a region 
much closer to the acceleration engine (~ 20-100 AU). The de- 
rived rotational motions appear to be in agreement with theoret- 
ical prediction s of models developed for magneto-centrif ugally 
launched jets dAnderson et all 120031 : iPesenti et aUl2004l) . sup- 
porting the idea that jets do transport angular momentum. 
These observations, however, only consider jets from evolved 
sources (Class II YSOs) in the last phases of their pre-main 
sequence evolution. 

If these signatures are evidence of jet rotation, then to 
establish the importance of this mechanism for star forma- 
tion we need to determine if such features are also seen at 
earlier evolutionary epochs, such as Class I YSOs. Current 
theories would predict this since inflow/o utflow is pre s ent as 
a disk form s arou n d the protos t ar (se e iPudritz et all [2007; 
IShang et all l2007h . iDavis et al.1 (120021) demonstrated that a 
small number of Class I YSOs harbour small scale H2 jets, so- 
called mole c ular h ydrogen emission line regions (MHEL - see 
IDavis et alJ (HH) extending some 2 - 3" from the central 
infrared source. These sources give us the opportunity to con- 
duct ground-based searches for and comparisons of the kine- 
matics of jet rotation between evolved T-Tauri and young Class 
I YSOs. 

Here we report on a pilot study, conducted on the VLT, 
which suggest that jets from Class I YSOs are rotating as they 
emerge from the central engine (§3), and present a simple kine- 
matic model which, at least for HH 26, supports this suggestion 
along with analysis assuming that the jet is launched magneto- 
centrifugally (§4). Finally, we comment on the relevance of 
these findings in the context of jet rotation studies (§5). 

2. Observations 

2.1. Source selection 

The jets from HH 26 and HH 72 were selected for this study. 
In selecting suitable Class I YSOs to search for jet rotation 
we are limited to using tracers longward of the optical. This 
is because sources this young are always found embedded 
within molecular clouds which also makes the use of adap- 
tive optics very difficult in the absence of a nearby optical 
guide star. The 2.12yum v=l-0 S(l) transition of H2is an ex- 
cellent tr acer of the dynamics of outflow s and jets in embedded 
sources (Ballv, Reipurt h & DavisLl2007l) . At present, the work 
of lDavis et al.ld200lLl2002l) represents the best source of YSOs 
with H2jets and two candidates from this work were selected 
to conduct our pilot study. HH 26 was chosen because it is rela- 
tively nearby and bright (in terms of available MHEL sources). 
Furthermore, it is a well studied object with published ancil- 
lary data to support our analyses. HH72 is less well studied 
and is much further away, but represents an unique opportunity 
to search for rotation in the jet from a more massive YSO. 

Our objectives with this exploratory study is to establish the 
presence or absence of evidence for jet rotation in these young 
sources before embarking on a more extensive study. 



observing log in Table[T]for details). 

ISAAC was configured using the SWS1-MR, medium res- 
olution spectroscopy mode and centred on the v=l-0 S(l) tran- 
sition of H2 at A — 2. 121 8 /mi . A spectral resolution of R=8900 
together with the 0.3" slit provided a velocity resolution of 
33.7 km s~' (or ~ 17kms~' per pixel). The instrument optics 
of 0.148" per pixel ensured good sampling of the emission 
structures (~ 2 — 3"). Each exposure through the long-slit 
spectrograph was 300 seconds. In some previous studies, slits 
were placed parallel to the flow and stepp ed across, building 
up a three-dim e nsion a l view of the j et (e.g . lDavis et all [2000; 
iBacciotti et ail 120021: IWoitas et all l2005h . However, the un- 
avoidable uneven illumination of the slits introduced spurious 
velocity shifts which had to be calculated and subtracted from 
each spectrum. In this work, the slit is placed perpendicular 
to the flow axis - at a position angle of 130° and 2" from 
the source of HH 26, and at a position angle of 0° and 1 .5 " 
from the star driving the HH 72 flow (see FigureQ]) - thus sam- 
pling the velocity field ac ross the flow in one exposure (see also 
Coff ey et all 120041 [2007 , who also observed in this manner). 

A standard ABBA observing sequence was adopted to re- 
move the sky background with the A and B telescope beam 
positions separated by 60" along the slit, sufficient to avoid 
any overlap of H2 emission and keep the source on the ar- 
ray at all times. The data were reduced using standard tech- 
niques and starlink software packages. Image straightening 
and wavelength calibrat ion was performed u sing OH sky lines 
(using the line list of [Rpusselot et al., 2000). The final wave- 
length calibration was accurate to ~ 0.1 A (rms dispersion 
of ~ 1.3kms _I ). The velocity scale was calculated relative to 
the vacuum wavele ngth of the H? v=l-0S( l) transition, given 
as 2.121833 /an by iBlack & van Dishoeckl (119871) . Velocities 
have been corrected for the velocity of the parent clouds^ and 
the Earth's relative motion on the nights observed. Position- 
velocity (PV) diagrams are constructed after aligning and com- 
bining the A and B beams. The seeing varied between 0.5" and 
1.0". 



3. Results 

3.1. PV diagrams 

In Figure[2]we present the position-velocity (PV) diagrams ex- 
tracted from the data for HH26 and for HH72. The slit offset 
is centred on the photocentre (brightest pixel) and is explic- 
itly assumed to define the outflow axis. The diagrams show 
that the jet velocities are quite high, of the order of -60 to 
-100 kms~' , while the velocity dispersion for HH26 is nar- 
rower (~ 45 km s" 1 ) than that for HH 72 (~ 80 km s" 1 ). 

The mean intensity-weighted radial velocity (measured on 
the axis) for HH72 is -56.3 + 0.3 kms 1 , co nsistent with 
previ ous H2 echelle spectroscopy for this object dDavis et all 
2001). Although observed at lower spatial resolution, so that 
the parti cular component p resented here was not properly re- 
solvecQ, Davis et all (l200lh did detect emission in their slit up 

1 -33 km s _I and +10 km s _I for HH72 and HH 26 respectively. 

2 The MHEL jet in bot h HH 26 and HH 7 2 were later identified us- 
ing Fabry-Perot imaging dDavis et"af , 2002). 
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Fig. 1. H 2 v=l-0 S(l) images of the HH 26 (left) and HH 72 (right) outflows (adapted from lDavis et al.Ll2002h . In each panel part 
of the larger scale flow is seen together with the YSO driving the outflow (from which the coordinates are referenced). The thin 
white boxes are representations of the slit placement and their position angles. The observations were targeted to the H2 jet close 
to the source with the slit placed perpendicular to the flow axis (see Section 2 for details). 



0.4 km s 1 . No pu blished H2 spectrosc opy towards the HH 26 
MHEL is known. IDavis et ail (|2000t) present echelle spec- 



to velocities of ~ -165 kms -1 . On source, they report a com- 
plex of emission features with components at +13, -40 and 
-129 kms 1 . As the slit for these data was placed ~ 1.5" from 
the source, we are probably detecting the same component re- 
sponsible for the -40 km s -1 peak. If so, then the jet either ac- 
celerates away from the source or was faster in the past. 

For HH26, the peak on-axis radial velocity is -87.5 + 
[2 spectroscc 
~aD (120001) 

troscopy of the H2 emission measured in the HH 26A and C 
outflow lobes (HH 26A is the lobe prominent in Figure Q] while 
HH 26C is further to the north-east and not shown here) with 
the slits placed along the outflow axis. Their data showed 
emission up to velocities of ~ -50 km s" 1 in HH26A and 
~ -90kms 1 forHH26C. 

For a slit placed perpendicular to the flow direction, a sig- 
nature of rotation is apparent as asymmetric contour lines in a 
PV diagram. Assuming axisymmetry, material moving towards 
and away from us at a given longitudinal radius is slightly more 
blue- and red-shifted with respect to material moving trans- 
versely to the line of sight. These differences are more evi- 
dent at large distan ces from the axis than closer in (e.g. see 



ICoffey et all 20041) beca use of projection and convolution ef- 
fects (IPesenti et alll2004T) . 

Asymmetry is present in both PV diagrams shown in Figure 
|2l the effect more pronounced in HH72 than in HH26. The 
figures also show that this feature is not only different in mag- 



nitude in the two objects but also of opposite sense (the direc- 
tion of the putative rotation determining which way the con- 
tour lines slope). Importantly, this indicates that it is unlikely 
that our results are due to an unforeseen or unaccounted for 
systematic error. 

3.2. Radial velocity measurements 

Figure|2]also shows the mean intensity-weighted radial veloci- 
ties measured for HH 26 and HH 72 along the slit and displayed 
as a function of distance from the axis. This velocity was cal- 
culated for each row along the slit by determining, 



<v> = 



2i h 



(i) 



where the summations are made over the whole spectrum. To 
determine the mean value and its uncertainty, each pixel in- 
tensity, /;, has its value randomly adjusted by either adding or 
subtracting its error, or not at all. This is repeated 1000 times 
for each pixel row thus establishing a distribution for (v). 

What is immediately striking in Figure [2] is the decrease 
in absolute velocity away from the axis of the flow. This 
has already been observed in a number of T-Tauri jets at 
angular resolutions high enough to resolve the flow width 
- of t he order of 100 AU at a similar distance from the 
source dBacciotti et alll2000tlWoitas et aUl2005l:ICoffev etaE 
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Fig. 2. PV and radial velocity diagrams for HH26 (top) and for HH72 (bottom) measured at a distance of 2" and 1.5" from the 
source, respectively. The left hand panels show the PV diagrams and the right hand panel shows the mean intensity-weighted 
radial velocity measurements (solid line); the dashed line shows the symmetry profile (the average of velocities either side of 
the axis) and approximates the profile in the absence of rotation. Offsets (in AU) are measured along the slit and relative to the 
pixel with highest flux. At the assumed distances to each object (see Table[TJ 1000 AU corresponds to 2.6" in HH26 and 0.7" in 
HH72. 



2004,|2007j). A decrease in poloidal velocity is expected from 
magneto-centrifugal launch theories, the Keplerian speeds in 
the launching zones of the disk wind naturally pr oviding a wide 
rang e of poloidal jet velocities across the beam dFerreira et all 
2006). In principle, it is possible that the outer portions of the 
flow are being slowed by interaction with the environment, 
however, we doubt this possibility given: (a) the jet is hyper- 
sonic and hence is not expected to en train much ambient ma- 
terial particularly close to the source dRav L 120001) . (b) there is 
no significant decrease in av erage poloidal velocity with dis- 
tance from the source (e.g. ICalzoletti et all |2008, who, for 
HH26, actually report an acceleration at the position of our 
observation), and (c) recent mass flux estimates along jets (e.g. 
Podio et all (2006) suggest very little if any entrainment. 

The dashed lines in Figure [2] are the symmetry profiles 
which would be obtained for unperturbed, non-rotating flows. 



It is computed from v'sym = (v+ + v_)/2, where v+ and v_ rep- 
resent the velocities measured at positions either side of and 
equidistant from the axis. 

The radial velocity data appear systematically displaced 
from the symmetry profile, with the direction of the displace- 
ment being different/opposite for the two objects. Such a dis- 
placement would be observed if the gas were rotating and, if 
so, the different directions indicate that the sense of rotation is 
different for the two objects. 

It is worth noting that this signature could also be gener- 
ated by asymmetric bow shocks, with the H2 streaming down 
the bow wings at different rates. However, within the present 
resolution limits we see no evidence o f bow shocks but ra ther 
a linear jet-like feature (see Fig. fTI and iDavis et all 12002)). If 
the H2 emission is tracing material shocked and entrained by 
the jet rather than the jet itself, then in the absence of alterna- 



follows we consider a simple model in which the observed PV 
diagrams are produced by a rotating jet. 



4. Modelling 

In this section we use analytical models designed to better un- 
derstand the relation between a rotating structure and the resul- 
tant PV diagram. Moreover, assuming that these structures are 
indeed rotating and using very general conservation laws for 
magneto-centrifugally launched jets, we derive some physical 
parameters from the observed features. 

4.1. PV diagram model 

Figure [3] is a schematic representation of an empirical model 
used to describe the section of a rotating YSO jet that emits the 
radiation collected by our slits. The jet section is represented as 
a geometrically thick disk (thickness - ze - zs) and is inclined 
towards our line of sight at some angle, 6. The jet section has 
an external radius Ro and is truncated internally at a radius Ri - 
inner axial holes are predicted in most models of X-winds an d 
disk winds (e.g. see iPudritz eTafl l2007t IShang et all I2007I) . 
The density distribution for the jet section is assumed to be a 
function of the radius, n(r), first increasing to a maximum and 
then decreasing smoothly outwards on a distance scale a. Such 
a distribution can be described by the function: 

n(r) oc r a exp(-r/a) (2) 

The material in the jet section possesses an axial velocity 
which we take to be the jet velocity, Vj e t(?"), and is rotating at a 
rate V^(r). The radial dependence of both these velocity com- 
ponents are expressed as power-laws, 

/ r Y P 

V") = v *>°(^j ( 3 ) 

VjetM = Vjet,o(^j (4) 

Note that, for convenience, the azimuthal velocity is nor- 
malised to the outer edge of the structure, while the jet velocity 
is determined relative to a central velocity defined at the inner 
truncation radius. With this simple formalism, a separate radial 
component v r (r) may also be added to reflect either expansion 
or contraction of the gas from or toward the axis. The velocity 
along the line of sight and at each point (x,y) of the plane (ie. 
edge-on) of the jet section is then given by, 

Kr) = v,(r)(^-v r (r)(2) (5) 

where x and y are coordinates on a Cartesian grid with the ori- 
gin coincident with r — 0. Allowing the structure to incline 
relative to the plane of the sky requires that the coordinate grid 
is transformed about the origin, 

x? — x 

y' = (ycos0) + (zsm8) (6) 
z' = (z cos 0) - (y sin ff) 




Fig. 3. Schematic of the structure used to model a rotating jet 
beam. The coordinate scheme, whose origin is centred on r — 
is also depicted. The x ordinate corresponds to the impact pa- 
rameter and z is the height or thickness of the object projected 
onto the plane of the sky. The Fb emitting portion of the jet has 
inner and outer radii truncated at Ri and Ro, respectively, and 
is rotating with a velocity v^. The transformation following an 
inclination of 9 about the x axis (towards the line-of-sight) is 
shown. 

Each velocity value at position (x, y, z) calculated by equation 
|5]is then assigned to position (x', y', z'). 

The velocity component in the inclined plane viewed along 
the line-of-sight thus becomes, 

v(r, 8) = v(r) cos 6 - Vjet(f) sin 6 (7) 

We calculate the H2 emission from the jet section entering 
our slit by assuming that the emitting gas is optically thin and 
dust free. The synthetic image thus produced is then dispersed 
in wavelength according to the instrumental characteristics of 
the spectrograph. In this way, we obtain a synthetic PV diagram 
which we can directly compare to our observations. 

4.1.1. Results 

Models whose jet velocity is constant across the jet do not pro- 
duce results which are representative of the data, regardless of 
the behaviour of the angular velocity component. Better rep- 
resentation is achieved if the poloidal jet velocity is allowed 
to decrease as a function of distance away from the axis. This 
alone shows that the jet velocities observed cannot be constant 
across the structure but must decrease with radius. 

However, the basic characteristic of a peak velocity profile 
offset from the peak of the symmetry profile is clearly present 
when rotation is introduced. Reversing the sense of direction of 
the rotation also displaces the velocity profile to the other side 
of the symmetry peak, indicating that it is the sense of rotation 
that is responsible for this characteristic. This compares quite 
nicely to what we see in Figure|2]and lends weight to our sug- 
gestion that the HH 26 and HH 72 jets are rotating in opposing 
directions. 

The shape of the PV diagram, and the corresponding peak 
velocity profile, is more sensitive to some model parameters 
than others. For a given inclination angle, models with steeper 
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Fig. 4. PV diagram models. The panels show PV diagrams along with the associated positions of intensity-weighted radial 
velocities (solid line) compared to a symmetric velocity profile (dashed line), i.e. with no rotation, beside them. The inclination 
of the flow from the plane of the sky is shown on the legends for each PV diagram. 



power laws for the jet velocity (parameterised by y in equation 
HI have the effect of introducing a prominent 'butterfly-wing' 
to the PV diagram. However, a weak dependence is seen with 
how the azimuthal velocity changes with radius, parameterised 
by f3, if y is fixed. In fact, the radial velocity diagrams remain 
remarkably similar for all values of /3 although they become 
more symmetrical as /3 becomes shallower. Any differences be- 
tween the models are subtle and most easily seen in the velocity 
diagrams. 

In Figure H] we show models as a function of inclination at 
fixed jet and azimuthal velocities (see Table [2] for the full list 
of parameters). The emission tracks across the PV diagram as 
the inclination increases, the jet velocity component becoming 
more prominent at higher inclinations. We also start to see limb 
brightened inner edges in the PV diagram, a consequence of the 
density structure (equation!!]), with the extent of the 'butterfly- 
wing' shape determined by the velocity power law (y). The 
radial velocity profile becomes more symmetrical at higher in- 
clinations. 

These models show just how instructive PV diagrams can 
be. Comparison with Figure|2]suggests that the jets are inclined 
to the line-of-sight at angles between 30 - 40°. 

No attempt has been made to actually fit a model to the data 
(the breadth of the parameter space precludes such a search). 
With the caveats already mentioned, the systematic changes 
seen in the intensity-weighted radial velocity diagrams are well 
represented by a jet that rotates about its axis. 

4.2. Analysis of the HH26 jet rotation signature 

In the following sections the results for HH 26 are analysed un- 
der the hypotheses that they do indeed indicate the presence 
of jet rotation, that the mechanism launching the H2 portion of 
the jet can be identified as a large scale magneto-h ydr odynamic 
(MHD) disk wind, (e.g. lKonigl & Pudritzl fcOOO) and lFerreiral 



Table 2. Parameters used in PV diagram model. 



Parameter 


Value 


Pixel scale 


0.148" pixel 


Distance 


400 pc 


Vertical Size^ 


1.0" 


Inner radius, Ri 


0.1" 


Outer radius, Ro 


2.5" 


Seeing (FWHM) 


0.5" 


Scale length (density), a 


0.5" 


Jet velocity, Vj et0 


150 kms- 1 


Azimuthal velocity, v^o 


3 km s _I 


Radial velocity, v r (r) 


0kms' 


Spectral resolution 


33.7 kms-' 


Exponent (density), a 


0.75 


Exponent (azimuthal velocity), {S 


0.5 


Exponent (jet velocity), 7 


0.1 



" Size of emission as it would appear on the sky. 



(2002)), and the quantities derived in our work can be used 
to constrain the properties of the jet acceleratio n region. The 
deriva tion follows closely the one described in IWoitas et al.1 
(2005) for the atomic jet from the classical T Tauri star RW 
Aur, seen at optical wavelengths. We do not conduct a similar 
analysis for HH 72 as there is insufficient information available 
in the literature to estimate values such as the mass accretion 
rate and mass flux. Furthermore, its greater distance from us 
than HH 26 makes the radial velocity measurements more dif- 
ficult to interpret. 

In the standard disk-wind model, the flow is assumed to 
be steady-state, axisymmetric and to satisfy the ideal MHD 
equations. In a cylindrical coordinate system (r,z,<p) with the 
star at the origin, the magnetic field lines identified by the 



plane at the so-called 'footpoint' (see Fig. 1 in IWoitas et alj 
(2005)). The disk and the anchored field lines rotate rapidly, 
and the centrifugal force launches fluid parcels from the disk 
surface out along the oblique open field lines. This in turn gen- 
erates a "magnetic torque" that brakes the disk, and extracts 
energy and angular momentum from it. In the acceleration pro- 
cess, the matter reaches a point (the so-called Alfven surface, 
at a height of a few AUs above the disk) where the inertia of the 
matter overcomes the magnetic forces, and the field is wrapped 
tightly generating a very important toroidal component. This in 
turn produces a magnetic force ("hoop stress") directed toward 
the axis, that collimates the flow. 

Using our radial velocity measurements, one can attempt to 
determine the so-called 'footpoint radius' of the wind, that is 
the location in the disk from where the observed portion of the 
wind originates, and the ratio between the poloidal and toroidal 
components of the magnetic field in the observed section of the 
jet. 

First, one has to compute the value of the poloidal and 
toroidal components of the velocity field, v p = v r e r + v : e z and 

respectively, at the observed location, that is at ~ 800 AU 
above the disk for HH26. These components can be derived 
from our rotation measurements simply as, 



Vp x v z = (v rad i + V rad2 )/(2cos0 

V<* = (v rad i - v rad2 )/(2 sin/), 



(8) 



where v ra <ii,2 are the measured radial velocities at two equidis- 
tant positions either side of the flow axis (with 1 referring to 
the northern part of the slit) and i = 90° - 6 is the inclination 
angle of the flow with respect to the line of sight. 

The derived values of these quantities are shown in Figure 
|5] The errors are propagated in quadrature from the measure- 
ment error of the radial velocity. We assumed a value of 40° 
for the inclination of HH 26, that follows from the combina- 
tion of tang ential proper motion measu rements (~ 70 km s _1 ) 
reported in Chrysostomou et al. (2000) for HH26 knot C, the 
closest to our observations, and our measured radial velocity 
(-86kms-'). 

The values for v p and are of course proportional to those 
shown in Figure [2] The poloidal velocity varies between about 
-75 and - 110 km s _1 while the toroidal velocity values span 
between 1 and 7 kms _I . One notes that while the poloidal ve- 
locity shows the expected behaviour, i.e., its magnitude gradu- 
ally decreasing with distance from the axis, the absolute value 
of V0 first increases reaching a plateau before decreasing again, 
while one would expect that it would have its maximum toward 
the jet axis. This is due to beam smearing and projection effects 
of the portion s of th e jet that are intercepted by the line of sight 
dPesenti et all 120041) . This effect may introduce an underesti- 
mate of the toroidal velocity by up to 20% in the two positions 
closest to the axis, which also affects the determination of the 
footpoint radius and of the magnetic pitch angle at the same 
position, as explained below. On the other hand one can see 
that the regular trend in v p and is interrupted in the last 2-3 
outermost positions, where the errors are largest. This indicates 
either an objective difficulty in determining the velocity, and/or 
the fact that the borders of the knot are interacting strongly with 
the surrounding environment, and the material emitting in these 
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Fig. 5. Results of the H2 rotation analysis for HH 26 assuming 
a flow inclination angle of 40° with respect to the line of sight 
and a central mass of 1 M . The values of the various quantities 
are obtained from the parameters of the wind observed at the 
distance from the axis defined by the x coordinate. From top 
to bottom: poloidal velocity; toroidal velocity derived from the 
radial velocity shifts; distance from the star of the wind foot- 
point, according to an ideal MHD disk wind scenario; absolute 
value of the magnetic pitch angle IB^/B^I at the observed loca- 
tions in the flow. 



positions is not participating to the general motion of the flow. 
For this reason the points at 1 - 1 .2" from the axis are excluded 
from the subsequent analysis. In conclusion, the most reliable 
estimates are obtained from data at positions 0.3 - 0.6" from 
the axis. 



4.2.1. Footpoint radius 

The third panel down in Figure [5] reports our estimate of the 
footpoint radius rn for the wind, i.e. the region in the disk from 
where the matter detected in our data at a distance r Q b s from the 
flow axis should originate. This can be derived from the com- 
bination of poloidal and toroidal velocities measured at r b s . 
In fact, in a magneto-centrifugal wind, the footpoint radius of 
the jet component located at r„b s from the axis with measured 
Vs and v p can be estim ated using a relationship provided by 
I Anderson et alJ d2003l) which is valid at large distances from 
the source: 



ro ,o. 7AU f^f 3 f v * -f 

UOAU/ UOkms- 1 / 



,/ Vp X-4/3/ M+ 

UOOkms- 1 / \1 



Mr. 



1/3 



(9) 



This simple equation is valid under the condition «: v p , 
which in the context of disk winds is achieved as soon as 
the Alfven surface is reached, at a few AUs above the disk. 
Thus at the sampled distances from the star, one can be confi- 
dent that this asymptotic regime has been reached by the outer 



tion of this luminosity is due to accretion one can infer that 
HH 26 has a mass of ~ 1 M Q . Once again, the errors are propa- 
gated quadratically from the measured quantities. 

Using equation d9) we obtain footpoint radii varying be- 
tween 0.5 to 6 AU from the protostar. Keeping in mind the 
caveat for measured velocities close to and far from the flow 
axis, one should take the extreme values with caution and ac- 
cept that a more reliable estimate of the footpoint radius for this 
jet lies between 2 and 4 AU for a star of 1 M (with an average 
error of 20%). 

Th ese values are consisten t with those of ICoffey et al.l 
(200% and IWoitas et al.l d2005l) when compared to the foot- 
points of jets emitting in optical lines. For those, which sample 
a portion of the flow nearer the axis, footpoint radii between 
0.5 and 2-3 AU were found. This is in agreement with the no- 
tion that the molecular wind should be produced in a region of 
the disk with low ionisation external to t he one from which op- 
tical jets orig i nate, a s found for DG Tau dBacciotti et all l2000: 
iTakami et aUl2004l) . 

It should be noted that the above calculation gives only the 
footpoint of the flow surface for which the toroidal velocity 
can be measured in the jet, and not the radius of the whole 
ejection region. In the disk-wind scenario, other colder com- 
ponents seen in molecular lines at longer wavelengths may be 
anchored to larger footpoint radii. 

4.2.2. Magnetic field properties 

In the hypothesis of disk- wind acceleration for the jet, our mea- 
surements allow one to easily find the ratio between the toroidal 
and poloidal components of the magnetic field, B^/Bp, in the 
jet at the observed location. This quantity indicates by how 
much the lines of force are wrapped on a given magnetic sur- 
face. The MHD models for jet acceleration prescribe that the 
collimation arises from the hoop stress generated by the in- 
crease of the toroidal comp onent of the magnetic field above 
the Alfven point (see, e.g. lKonigl & Pudritzl d2000l) '). It is thus 
interesting to check if the observations give any indication for 
the magnetic field configuration consistent with this property in 
the region of the flow just above the collimation zone. The ratio 
B$l B p is foun d using another conservation law of general dis k- 
wind theory dKonigl & Pudritzl 120001: [Anderson et afll2003l) : 

(vv-^pJ/r = Q , (10) 

where Qo is the disk angular velocity at the footpoint. 

Introducing our measured velocities at r — r \, s and using 
the corresponding ro to calculate Qo for a Keplerian disk, we 
find BtyjBp for any r = r \, s . The results are shown in the bottom 
panel of Figure|5]in absolute terms (the sign of the ratio is con- 
sistent with that of the velocity components and the assumed 
geometry). Once again caution should be taken for points both 
close and far from the axis. Good estimates are probably given 
by the three positions between 0.3 - 0.6", for which we find 
\B</,/B p \ varying from 10 to 20. 

Thus, the toroidal component of the magnetic field appears 
to be dominant, which supports the suggestion that protostel- 
lar outflows are collimated magnetically, as prescribed by the 
MHD acceleration models. 



Following magneto-hydrodynamical models, the angular mo- 
mentum transported along each flow surface of the jet is a 
constant along the surface and has both a kinetic and mag- 
netic component that always add in sign. The total angular mo- 
mentum transported by the flow can then be calculated by in- 
tegrating the contributions of all the flow surfaces across the 
jet. For the two systems for which this has already been done 
(namely, DG Tau and RW Aur), the angular momentum carried 
away was estimated to be between 60% and 100% of the angu- 
lar momentum that the inner disk has to lose to accrete at the 
observed rate dBacciotti et all 120021; IWoitas et al.l 120051) . The 
fundamental implication of this is that jets are likely to be the 
major agent for extracting excess angular momentum from the 
star/disk system which allows matter to flow through the disk 
and accrete onto the star. 

We cannot undertake a similarly accurate calculation for 
HH 26 as the internal distribution of the molecular gas in the 
jet, and therefore the location of the most internal flow surface 
of the H2 -emitting gas, is not precisely known. This affects 
the calculation of the integral of the angular momentum and in 
particular the deter mination of the magnetic contribution (see 
IWoitas et all 120051 . Nevertheless, it is possible to approximate 
the kinetic contribution, a lower limit to the actual angular mo- 
mentum transported, and compare this part with the angular 
momentum that has to be lost by the disk. The kinetic angular 
momentum can be expressed as, 

ijet = J"(^)pv P ndS « (rv^)M jet , (11) 

where S is the surface of the observed section of the jet, p is 
the total mass density, (rv$) is the average value of the specific 
angular momentum in the observed region and Mj et is the jet 
mass flux in H2 . 

For HH 26, the H2 mass flux r ate in the internal knot co nsid- 
ered here has been computed in Ant oniucci et al.l J2008) who 
derive a value between 2 - 5 x 10~ 8 M Q yr _1 . The same authors 
also estimate a mass accretion rate of 8.5 x 10 -7 M yr~' from 
the luminosity of the Bry line. Assuming that the toroidal ve- 
locity measured at r D b s = 0.44" , = 2.5kms~' is representa- 
tive of the average over the H2 emitting region, one obtains for 
a distance of 400 pc, L jet ~ 2 x 10~ 5 M yr _1 AU km s _1 . For 
comparison the total angular momentum transported by the op- 
tical jet in the RW Aur system was estimated to be ~ 2.5 x 10~ 5 
M Q yr~' AU km s _1 in each lobe. 

It is now interesting to compare this with an estimate of the 
angular momentum the disk has to lose for the star to accrete at 
the observed rate. The excess dis k angular momentu m can be 
found through the expression (cf. Woit as et alll2005l) : 

^disk = (M^cc + -Mflow) ni.ext ^K.ext ~ -^acc ''O.in ^'K,in (12) 

where M acc is the mass accretion rate onto the star, Mfl ow is 
the sum of the mass loss rates in the two lobes of the whole 
flow (atomic and molecular components), r^in.ext are the inner- 
and outermost footpoint radii of the whole flow and VK,in,ext 
are the corresponding Keplerian velocities. Both magneto- 
hydrodynamic models and observations prescribe that Mfl ow = 
0.1 M acc , and that r ,i n ~ 0.03 AU from the star. 



^disk ~ 6 x 10 M G yr~ J AU km s . Assuming that HH 26 has 
a counterjet that possesses similar physical properties, we find 
that the kinetic angular momentum transported by the bipolar 
H2 jet already amounts to about 70% of the angular momentum 
that has to be extracted from the disk. Considering that the ki- 
netic component of the molecular jet gives just a lower limit to 
the actual total angular momentum in the flow, our result lends 
further support to the hypothesis that jets are the main agent for 
extraction of excess angular momentum in the disk. 

4.3. Caveats and alternative explanations 

The results and analysis presented here seem to indicate that, 
for HH26 at least, the observed portion of the jet originates 
from an extended region of the circumstellar disk. We cannot 
exclude, however, t h e pres ence of an inner X-win d or stellar 
wind (cf. IShu et alj (l2000h and ISautv et alj ( 120031) "). as at the 
resolution of our observations we only probe the external flow 
surfaces. In addition, it is not yet proven that flow lines can 
be entirely traced back from the observed jet locations to the 
footpoints in a real jet. In fact, once the region beyond the ac- 
celeration zone is reached, various magneto-fluid instabilities 
can complicate the geometry of the field lines and a direct con- 
nection to the footpoints may be lost. 

Furthermore, a number of recent studies have proposed al- 
ternate explanations for the observed velocity asymmetries. In 
fact, there are aspects that still await clarification, the most in- 
triguing of which to date is the recent study of the kinematics 
of the RW Aur disk that suggest a rot ation sense for the disk 
opposite to the one of the bipolar jet (Cabrit et al., 2006). At 
the moment, therefore, alternative explanations for the obser- 
vational evidence are actively searched for. According to re- 
cent studies, velocity asymmetri es may als o be produced by 
interaction with a warped disk dSokerL 2005h. or in asymmet - 
ric shocks generated by jet precession (Cerqueira et al., 2006). 
However, in the ir simulations of rotating and/or precessing jets, 
ISmith & Rosenl (120071) find no evidence of rotation signatures 
being mimicked by a precessing jet, although they do warn that 
any rotation that is present will quickly be dissipated as the jet 
expands. Interestingly, they predict that jet rotation introduces 
instabilities in the flow which develop into a "swarm of bow 
features", rather like those observed ahead of our slit place- 
ments (see Fig. [TJ. Such interpretations need to be able to ex- 
plain the statistics (in practice 100% of examined cases) with 
which velocity asymmetries are found, and within the appro- 
priate range predicted for rotation. 



5. Conclusions 

We present observations of the H2 emission from jets of two 
embedded Class I sources, HH 26 and HH 72, and search for 
velocity asymmetries in the spectral lines from across the jets 
that are compatible with rotation. Our results are consistent 
with those signatures found at optical wavelengths for T Tauri 
jets, but for objects at an earlier evolutionary phase. 

The position-velocity diagrams for HH 26 and HH 72 show 
a tilt in the intensity distribution of the PV diagrams, i.e. an 
asymmetry of the jet velocities with respect to the axis which 
can be interpreted as rotation of the gas. The observed tilt is 



PV diagrams using an analytic model for the H2 radiation emit- 
ted by a cylindrical jet observed by a spectrograph slit. The ob- 
served skew in the PV diagrams can be reproduced assuming 
that the jet material possesses a toroidal velocity of the order of 
a few kms~' which decreases as a simple power-law with the 
jet radius, together with a jet velocity that also decreases with 
distance from the jet axis. 

Under the assumption that the flow is indeed rotating, and 
that it is launched magneto-centrifugally, one can apply very 
general conservation equations valid for such jet acceleration 
mechanisms to retrieve information on the physics of the jet it- 
self. One can show that the part of the flow under consideration 
may originate in an annulus of the disk located at a distance 
from the star of 2-4 AU for HH26. This agrees with the no- 
tion that the molecular wind should be produced in a region of 
the d isk external to the one from which the optical jets origi- 
nate dTakami et al.L 120041) — up to 2-3 AU in the disk accord- 
ing to previous rotation studies. In a similar way it is possible 
to show that the toroidal component of the magnetic field in 
the observed portion of the flow is dominant over the poloidal 
component, in analogy with what has been found for optical 
jets and in support of a magnetic collimation mechanism. 

On the basis of our measurements and previous determi- 
nation of mass loss rates, we estimate that the kinetic an- 
gular momentum transported by the HH26 jet is about 2 x 
10 _5 M o yr _1 AUkms~' . This is a lower limit to the total an- 
gular momentum transported, which also includes a magnetic 
contribution (not derivable from our measurements), but al- 
ready amounts to 70% of the angular momentum that has to 
be lost by the disk for the star to accrete at the observed rate. 
A similar estimate for HH72 has not been possible, due to un- 
certainties introduced by its large distance and the lack of an 
estimate of the accretion rate onto the source. 

These results are in agreement with most model prescrip- 
tions of protostar formation according to which the jets are 
produced in the system as soon as the disk begins to accrete 
material onto the central protostar. This occurs because of the 
need to remove excess angular momentum from the central ac- 
creting system. Further observational studies, using different 
tracers and conducted on a larger number of jets as well as 
their associated disks, need to be made towards objects at vari- 
ous evolutionary phases in order to prove whether rotation is a 
general property of stellar outflows. 
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